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ABSTRACT 

A  Carousel  Remotely  Piloted  Vehicle  (RPV)  Recovery 
System  consists  of  a  vertical  mast  which  supports  a  horizon- 
tal member  that  is  free  to  rotate.   The  approaching  RPV  is 
caught  by  vertical  cables  suspended  from  the  horizontal 
member  and  the  kinetic  energy  of  the  RPV  is  dissipated 
through  the  motion  of  the  recovery  system.   This  thesis 
presents  a  simplified  dynamic  analysis  to  describe  the 
motion  of  an  RPV  after  impact  with  a  recovery  system  simu- 
lating the  carousel  system.   The  equations  of  motion  for 
the  RPV  were  obtained  from  Lagrange's  equations  and  modeled 
using  the  Continuous  System  Modeling  Program  on  the  IBM 
360-67  computer  at  the  W.R.  Church  Computer  Center,  NPS . 
The  results  showed  that  damping  was  required  on  motion  of 
both  the  horizontal  member  and  the  suspended  cables  in  order 
to  prevent  possible  damage  to  the  150  pound  RPV  at  an  as- 
sumed 50  knot  approach  speed.   A  sensitivity  analysis  was 
performed  by  varying  the  system  design  parameters  and  the 
results  are  presented  in  the  body  of  the  thesis. 
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LIST  OF  SYMBOLS 

D         -  Rayleigh  dissipation  factor 

ft        -  Feet 

g         -  Acceleration  of  gravity,  (32.2  ft/sec) 

I         -  Moment  of  Inertia  of  horizontal  member  about 
o 

its  pivot  point,  (it-lb-sec1-) 
K         -  Spring  constant  of  torsional  spring  located  at 

pivot  of  horizontal  member  (ft-lb) 
lb        -  Pounds  force 
M         -  Mass  of  RPV,  (ft-lb-sec2) 

q         -  Generalized  coordinate 

rad       -  Radians,  angular  measure 

R-j_        -  Length  of  horizontal  member,  (ft) 

R2        -  Length  of  massless  recovery  rod  corresponding 

to  RPV  distance  from  horizontal  member  following 
recovery  (ft) 

sec       -  Second,  unit  of  time 

T         -  Kinetic  Energy,  (it -lbs) 

V         -  Potential  Energy,  (ft-lbs) 

V-Rpv      -  Velocity  expression  for  RPV  (ft/sec) 

x,y,z,     -  Fixed  axis  coordinates 

x',71 ,z'   -  Moving  axis  coordinate 

^  -   Dimensionless  damping  ratio 

0  ,  (p  t  iff     -   Coordinate  system  angles 

0max     -  Maximum  magnitude  of  angle  0 
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(p         -  Initial  condition  on  angular  velocity 
(jq  -   Undamped  natural  frequency,  (l/sec) 

-  Column  vector 

-  Matrix 

-  Determinant 

Note:   Dot  over  a  symbol  represents  differentiation  with 
respect  to  time;  i.e.,  x=dx/dt 
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I.   INTRODUCTION 

Extensive  efforts  on  Remotely  Piloted  Vehicles  (RPV) 
have  been  underway  for  many  years,  covering  a  wide  range  of 
applications  and  missions,  Ref.  1.   It  is  only  recently, 
however,  that  both  the  need  and  the  technological  capability 
for  mini-RPVs  have  coalesced  to  the  point  that  the  develop- 
ment of  these  systems  has  become  imperative.   A  mini-RPV 
is  generally  defined  as  an  RPV  that  weighs  less  than  200 
lbs  mission  loaded. 

An  essential  part  of  an  RPV  system  is  the  recovery 
system.   The  Army  has  recently  experienced  several  mini- 
RPV  losses  due  either  directly  or  indirectly  -to  the "recovery 
method.   The  methods  used  to  recover  mini-RPVs  vary  from 
catching  the  RPV  with  an  arresting  cable  similar  in  prin- 
ciple to  that  used  on  an  aircraft  carrier,  to  guiding  the 
mini-RPV  directly  into  a  vertical  net  which  is  designed  to 
flex  and  absorb  the  kinetic  energy  of  the  RPV.   A  completely 
satisfactory  recovery  system  has  not  yet  been  established, 
hence  studies  into  alternate  recovery  systems  are  of  con- 
tinuing interest  to  the  Army. 

The  purpose  of  the  analyses  described  herein  was  to 
model  the  dynamics  of  one  alternate  recovery  system  candi- 
date, the  "carousel"  method.   The  carousel  method  was 
initially  described  as  a  candidate  for  a  recovery  procedure 
in  a  concept  study  performed  by  Teledyne  Ryan  for  the  Army, 
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Ref.  2,  but  the  initial  concept  formulation  did  not  include 
significant  dynamic  considerations.   Therefore,  the  analyses 
conducted  in  this  thesis,  which  are  primarily  oriented  to- 
wards system  dynamic  response  behavior,  should  serve  in  the 
role  of  providing  additional  information  for  concept  feasi- 
bility considerations. 

The  carousel  recovery  system,  shown  in  Figure  1,  was 
conceptualized  in  this  analysis  as  being  composed  of  a 
vertical  mast  that  supported  a  horizontal  member  which  was 
allowed  to  rotate  about  the  vertical  mast.   Hanging  from 
the  horizontal  member  are  cables  that  will  be  caught  by 
the  wing  tips  of  the  approaching  RPV.   The  RPV's  kinetic 
energy  will  be  dissipated  by  the  recovery  system  through 
various  means  that  will  be  evaluated  in  this  analysis. 

The  RPV  must  be  brought  to  rest  without  allowing 
excessive  RPV  swing  since  excessive  swinging  motion  could 
result  in  the  RPV  hitting  the  recovery  system  structure. 
This  consideration  was  of  paramount  concern  when  evaluating 
the  merits  of  recovery  system  parameter  variations. 

For  the  purpose  of  this  analysis  the  mini-RPV  to  be 
recovered  was  assumed  to  be  the  AQUILA,  which  is  a  U.S. 
Army  mini-RPV  developed  by  Lockheed  Missiles  and  Space 
Company,  Ref.  1,  and  is  sketched  in  Figure  2.   It  has  the 
following  nominal  characteristics:   wing  span,  12  ft; 
weight  150  lb;  approach  speed  50  knots;  maximum  loading, 
6g  vertical  and  axial,  3g  lateral. 
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II.   ANALYSIS 

A.   ASSUMPTIONS 

Simplifying  assumptions,  as  defined  below,  were  made 
when  establishing  the  equations  of  motion  using  the  pro- 
cedures suggested  "ay   Lagrange's  equations,  Ref.  3>  in  order 
to  establish  the  basic  principles  of  the  carousel  recovery 
system  dynamics.   Inclusion  of  system  details  such  as  re- 
covery boom  elasticity  and  RPV  mass  mome'nts  of  inertia  were 
considered  as  items  which  might  tend  to  mask  the  basic 
question  of  feasibility  at  the  present  stage  of  analysis. 
Under  these  assumptions,  the  results  supported  the  viability 
of  the  carousel  recovery  method  and  a  basis  was  established 
upon  which  further  analyses  may  be  performed  when  the  sim- 
plifying assumptions  are  removed. 

The  recovery  system  was  initially  assumed  to  be  con- 
servative, that  is,  friction  losses  were  ignored  and  as 
soon  as  the  RPV  contacted  the  recovery  system  the  engine 
was  secured.   This  assumption  was  subsequently  changed  by 
the  addition  of  damping  to  the  recoverjr  system  configuration, 
The  horizontal  member  was  assumed  to  be  rigid.   A  careful 
choice  of  material  and  structural  characteristics  for  the 
horizontal  member  can  minimize  its  dynamics  and  allow  it  to 
be  approximated  by  a  rigid  member.   The  cables  hanging  from 
the  horizontal  member  were  replaced  by  massless  thin  re- 
covery rods.   The  RPV  was  treated  as  a  ooint  mass  and  was 
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captured  by  one  recovery  rod.   Treating  the  RPV  as  a  point 
mass  simplified  the  coordinate  system  used  in  the  analysis 
by  allowing  the  mass  moment  of  inertia  of  the  RPV  about  its 
principal  axes  to  be  neglected.   The  RPV  was  also  assumed  to 
approach  perpendicular  to  the  plane  formed  by  the  horizon- 
tal member  and  the  vertical  recovery  rod  that  caught  it. 
The  approach  velocity  of  the  RPV  was  assumed  to  transfer 
undiminished  to  the  recovery  rod  giving  the  rod  an  initial 
angular  velocity  only  in  the  0  direction.   Reference  should 
be  made  to  Figure  3  where  the  orientation  of  the  generalized 
coordinates  0,  0  and  ij/   is  portrayed.   These  coordinates 
are  a  form  of  Euler  angles,  independent  of  each  other,  that 
together  completely  specify  the  position  of  the  RPV  relative 
to  the  fixed  reference  frame  x-y-s .   First  0  defines  the 
rotation  of  the  moving  coordinate  system  x'-y'-z1  about  the 
fixed  axis  z.   It  describes  the  rotation  of  the  horizontal 
member  in  the  x-y  plane.   Next  0  defines  the  motion  of  the 
RPV  about  xf  within  the  moving  coordinate  system.   This  is 
portrayed  in  Figure  3  as  the  motion  of  the  RPV  projected  on 
the  y'-z1  plane.   Finally,  "J/^  depicts  the  remaining  motion 
of  the  RPV  after  it  rotates  through   0  and  0  . 

The  system  identified  in  Figure  3  is  the  reference  or 
base  configuration  that  was  assumed  for  the  analysis  based 
on  judgement  of  the  RPV's  weight  and  size.   In  this  refer- 
ence configuration  the  horizontal  member,  B.    ,  was  23  ft 

long  and  made  of  6-inch  aluminum  I-beam  having  a  mass 

o 

moment    of   inertia   of   1159   ft-lb-sec      about   -pivot    "0". 
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Figure   3      Simplified  Recovery  System 
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The  recovery  rod  R_,  was  6  ft  long.   The  RPV  of  mass  4.66 
slugs  (150  lbs)  was  suspended  at  the  lower  end  of  R~.   With 
an  assumed  approach  speed  of  50  knots  (84.4  ft/sec),  the 
initial  angular  velocity,  0  ,  of  the  massless  rod  Rp  was 
14  rad/sec  according  to  the  relation  of  equation  (l) : 

\          Ap-oroach  velocity  (ft/sec) 
$0  =    Hg  (fb)     (1) 

A  linear,  torsional  spring,  of  spring  constant  K,  lo- 
cated at  "0",  the  pivot  of  the  horizontal  member,  was 
initially  used  to  resist  the  rotation  of  the  horizontal 
member.   without  any  damping  in  this  reference  configura- 
tion it  was  further  assumed  that  the  spring's  stored  energy 
could  be  eliminated  from  the  system  through  a  ratchet  re- 
lease mechanism,  after  the  horizontal  member  had  swung 
through  its  initial  maximum  displacement.   The  ratchet 
mechanism  was  not  modeled,  as  it  was  felt,  and  later  con- 
firmed by  the  results,  that  the  critical  motion  of  the  RPV 
would  occur  shortly  after  impact  with  the  recovery  system. 
Simulation  of  the  ratchet  mechanism  would  have  no  effect 
on  these  results. 

E.   OVERVIEW  OF  ANALYSIS 

The  general  approach  of  the  analysis  was  to  model  the 
RPV  motion  after  engagement  with  the  recovery  system  using 
the  generalized  coordinate  system  @-0- 1//"  indicated  in 
Pig.  3.   If  the  assumptions  of  a  rigid  horizontal  member 
and  a  point  mass  RPV  representation  were  removed,  a  minimum 
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of  five  more  generalised  coordinates  would  have  been  re- 
quired.  The  equations  of  motion  for  the  three  degree  of 
freedom  system  were  obtained  using  Lagrange's  equations 
for  a  conservative  system.   From  these  equations  the 
acceleration,  velocity  and  position  of  the  RPV  as  a  function 
of  time  were  obtained  in  component  form  along  each  coordi- 
nate direction  through  the  use  of  the  Continuous  System 
Modeling  Program  (CSMP) .   The  input  parameters  were  then 
varied  and  their  effect  on  resulting  dynamic  motion  studied. 

Changes  to  the  basic  recovery  system  which  included  the 
addition  of  damping  terms  were  later  entertained  as  a  con- 
sequence of  initial  results. 

C.   DETAILED  ANALYSIS 

1.   Lagrange's  Equation 

The  equations  of  motion  were  established  using 
Lagrange's  equation,  Ref.  3.   Three  equations  of  motion 
were  obtained  for  the  system,  since  there  were  by  assump- 
tion only  three  generalized  coordinates:   the  coordinates 
0,  0    and  y   ,  as  defined  in  Fig.  3.   Lagrange's  equation 
in  the  conservative  form,  without  any  generalized  forces,  is 


_d  [  _^_T  \   _^_T    _Vv 

at  \  w,  r  aou  +  ^  (2) 


+  -£-1  =  0;      i  =  1,  2,  3, 
•i  /   --*i    d  -i 


whe  re 


ci-j   =   generalized  coordinate  (q-,=  0  ,Cp=0  ,q^=T^*) 
a  •   =   time  derivative  of  q.  coordinate 
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T   =  TUi-j^c^)  =  Kinetic  energy  (ft-lb) 
V   =  V  (q.)  =  Potential  energy  (ft-lb) 

The  kinetic  energy  was  defined  relative  to  the  inertial 
reference  frame  by  the  following  expression: 

t  =  *io  e  2+iMEPV2  (3) 


wiie  re 

I     =  Moment  of  Inertia  of  the  horizontal  member 
o 

about  pivot  point  "0"  (ft-lb-sec^) 

I.I     =  Mass  of  RPV  (slugs) 

0     =  Angular  velocity  in  0  coordinate  direction 

(r ad/sec) 
Y-c-n-r     =  Velocity  expression  for  RPV  derived  in  the 

next  section  (ft/sec) 

The  potential  energy  included  the  gravitational  term  acting 
upon  the  RPV  mass  and  the  stored  energy  in  the  torsional 
constraint  spring.   Equation  4  expresses  the  potential 
energy  relative  to  the  inertial  frame  of  reference  as 
follows: 

V   =  MgR2(l-cos0cosT/r)+iK02  (4) 

where 

Rp   =   length  of  recovery  rod  corresponding  to  RPV 

distance  from  the  horizontal  member  following 
recovery,  (ft) 
K  =   spring  constant  of  spring  at  pivot  of  horizontal 

member  (ft-lb) 
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Since  its  exactness  is  critical  to  the  evaluation  of 
Lagrange's  equations  of  notion,  the  RPV  velocity  is  treated 
separately  in  the  ensuing  section  using  the  generalized 

coordinates. 

2.   Derivation  of  RPV  Velocity 

To  derive  an  expression  for  the  RPV's  velocity 
after  impact  with  the  recovery  system,  the  position  of  the 
RPV  within  the  moving  x'-y'-z'  coordinate  system  as  a  func- 
tion of  the  generalized  coordinates  0,  iff  was  transformed 
to  the  fixed  inertial  reference  system  x-y-z.   The  time 
derivative  of  the  RPV's  position,  expressed  as  a  function 
of  the  three  generalized  coordinates  0,0  and  ljf   and  ref- 
erenced to  the  fixed  inertial  system,  represented  the  ab- 
•solute' velocity  of  the  RPV. 

The  RPV's  position  within  the  moving  frame  x'-y'-z' 
was  expressed  as  follows: 


t      _ 


y    = 


_r0    sin  iff 


&2    sin0cos7//" 
-R2    cos0cos  llf 


(5a) 

(5b) 
(5c) 


Transformation  to  the  fixed  inertial  coordinate 
system  may  be  stated  as: 


r    ~\ 

X 

< 

>  = 

y 

L  j 

z       = 


cos0  -   sinQ 
sin0       cos  0 


r 


-R   s±nlpr 


*> 


r 

R   cos0 


+  < 


Rpsin0cos"tyH         R1sin0 


v_ 


j 


z'  -R9cos0  cosier 


(6a) 
(6b) 
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In  an  alternate  form,  the  three  components  expressing  the 
absolute  position  of  the  RPV  mass  following  recovery  en- 
gagement may  he  expressed  as: 


x     =     R-,  cosQ  -R0sinl//r'cos  0 -RpSin0  cosl//"sin0 

y      =     R   cosQ  -R   sinl/TsinQ  +R   sin0  cosl//~cos  Q 
z      =      -R   cos0cos  Ur 


(7a) 

(7b) 
(7c) 


The  absolute  velocity  of  the  RPV  was  obtained  by  taking  the 
first  time  derivative  of  equation  7,  thereby  yielding: 


x     = 


y     = 


f-R1sin@  +R2(sinl//"sin0  -sin0cosl/Tcos@  )1  Q 
-(R2cos0  cosl/^sinQ  )0  +R9(  sin0  sinl/rsin0  - 
cos^cosQ  )  ijf  (8a) 

R..COS0  -R2(sinl/^cos0  +sin0  cosl/rsin0  )"]  Q 

+  (R-.cos0cos'l/rcos0  )0  -R0( cosl/fslnQ  +sin0sinl/r 

cos0   )tjf  (8b) 

z      s     R9  [(sin0cosl/r)0 +(cos0sinl/r)   l/r]  (8c) 

The  expression  desired  for  use  in  Lagrange's 

2 

equation,  VRpV  ,  was  obtained  by  squaring  the  components 

separately  then  adding  them;  i.e., 


V 


xiPV 


2  .2  .2 

:  +y  +z 


(9) 


•  2   -2      -2 
The  expressions  for  x  ,  y   and  z   are  equations  10a,  b 

and  c . 
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2  A2 


F-R-jSinQ  +R0(sinl/Tsin  0-sin0cos"|/rcos  0  )~|   "  q 


2_    2 


2^:2 


+R2    cos"0cos^^sin'~0  0 

9  9*9 

+r2   (sin0sinl/rsin0 -cos  *^cos0  )    l/r 

-2R2COS0cosVrsin6  |-R-i  sin@+R?(sin^rsin@  -sin0 

cosl/rcos0  )]  g  0 

+  2R   (sin0sin7/rsin0  -cosl/Tcos@  )  ["-R-,  sing  +R    ( sin  iff 

sin0 -sin0cosTl/rcos0  )  0  Ur 

-2Rp   cos0cosT/TsinQ  (sin0  sinT^sing  -coslJfcosQ  )0t/t 

(10a) 


.2 


v 
1/ 


9     '2   „    2 


JR,  cos0  -R9(sinl/Tcos0  +  sin0  cos^sin  0  jl-      0  +R 


L2 


2  2  2*2 

cos    0cos"~l/^cos   0  0" 


9     ,  9 


+R2     rcosl^rsin0+sin0sinT//r,cos0"|  "  }f/~ 

+2RpCOS0cosl/rcos  0  [r.,  cos0  -R2(sin"^Tcos  0+sin0 
cosl/Tsin0  )]0  0 


-2R 


cos"^Tsin0  +sin0sin"j^rcos  0]    |R,  cos0  -R„( sinl/r 

cos0  +sin0cosT^*sin0    0"|/r 

2  r 

-2Ro'"cos0cos  ijfcos  Q    cosT^*sin0  +sin0sin-|/r 


cos 


(10b) 
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9  O    r  O  9  9  9  O  '  O 

z        =     R_     bin  0  cos   iff  0   +  cos  0  sin   '^T^lf  +2sin0cos0 

sinl/^cosl/T  0  T/rl  (10c) 

•  2   >2   •  2 

When  x  ,  y  ,  z4"  were  combined,  several  cancellations 

2 
occurred  and  the  final  expression  for  V^,™,   simplified  to: 

VRpv2  =  [\2+R22   sin2T/r+R22sin20  cos^^R-^sinl/r]  Q2 

+R2  cos2Yr02+R22i(^2+2R2cos0cosi/r(R;L-H2sinl/r)  0  0 

2       ' 
+  (-.2R1R2sin0sin"r +2R2  sin0)6^ 

(11) 

3 .   Derivation  of  Equations  of  notion 

a.  Q    Equation  of  Motion 

2 

Now  that  an  expression  for  VRPy   has  been 

obtained,  Lagrange's  equation  can  be  evaluated  according  to 
section  II.C.l.   The  equation  was  evaluated  for  each  gener- 
alized coordinate  0  ,  0  and  ijf   .   The   0  equation  of  motion 
was  obtained  by  evaluating  the  Lagrange's  equation  using 
the  0  generalized  coordinate. 

Again: 

9  9 

Kinetic    energy,    T=J-I0  0   +iWRpy  (!3) 


Potential   energy,    V=MgRp  Il-cos0  cosT/n  +■§£  0 


(14) 


1  =  °  (15) 
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}Q  -  Ly  (16) 


d/9T 

cR; 


l^]=I0G  +M  [(R12+R22sin2-^+R22sin20cos2T^-2R1R2sinYr)0 
+R2cos0cosl/^(R-L-R2   sinl/r)0  +  (-^-^23^03^1^+ 


2 


R2    sm 


(pit 


2  2  2 

+  (2R2    sinl/rcos"l^i^+2R2    sin0cos0cos   7/^0 

-2R   2sin  0  cosljf  sinljfljr^R-^R 2c os -[{f^lr)  Q 

'  2 
-R2sin0cosl/r(R1-R2sinl/r)  0   -R2sinl/rcos  0(R;,-R 


sinl/O  0-j^r 

2  '    •  •     • 

-R2   cosT/^cos0cos7//"l/r0-RiR2Cos0  sini£r0  Yr-R]_R2 

•  2 
sin0cosl/r-|^r 

2  ' 

+R2   cos0  0i^r  /17x 


Substituting  equations  15,  16  and  17  into  equation  12  and 

2 

dividing  through  by  (MR~  )  yields  the  following  ^equation 

of  motion: 
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ID 

+  (Tp-sinT/T)2+cos2T/rsin20  )  Q 


mr22    xt2 

Rl 
+  cos0  (^-cosi^/--^sin2  7/r)  0 

2^-1  .     . 

+  sin27//*cos  0  -2—cosT/O  i/r  0 


2  '     • 

+  sin20cos   lJ/~  0  Q 


+  CO30  (-2^-sinl/T-cos2T/r+l)  0  ^ 

R 

1  •  2 

-sin0(TT-cosi^--5-sin2T/r)  0 


Rl  '  2 

-  ^sin0cosl/ri/^ 


+  ^2      =0 


lotto 
c. 


(18) 


b,  0  Equation  of  Motion 

Similarly  the  0  equation  of  motion  was  ob- 
tained from  the  following  Lagrange's  equation: 


a  (25  \       2r    +    av  =  0  (19) 

30/  ^0  <)0 


where 
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dtUp 


R2   cos  T/T0  +R9cos0  cosT/r(R1-R2sinT/T)  0 


-2R22cosl/Tsin'l/r-^0-R2sin0cosl/^(R1-R2sin-|/r)0@ 

-R2sinT/rcos0  (R1-R2sinl/r)@-^r  -R22cos0  cos2^]/^] 

(20)" 


1£  =  r,l[R22sin0cos0cos2Yr02-R2sin0  cosl/rCR^^-RpSin-^r)^^ 


30 


+  (-R-jR2Cos0  sin"l/r+R2   cos0)6l/r] 


(21) 


3V 


90      2 


(22) 


Substituting  equations  20,  .21  and  22  into  equation  19  and 

2 
dividing  through  by  (MR-  )  yields  the  following  0  equation 

of  motion: 


R- 


^— cos0  cosier -v;cos0  sin2l/r 


e 


+  C03 


V0 


-sin2l/^  1/^0 
-cos0(cos2T//-+1)0t/t 
--|sin20cos2"l^  Q2 


+  g-sin0cosl/r=   0 


(23) 
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c.  iff   Equation  of  Motion 

The  Tj/~  equation  of  motion  was  similarly  obtained 
using  the  following  Lagrange's  equation: 


a/  9t  \  _  3t      3v  _  0 

^l  j^i    #-   a^~  (24) 


where 


al(j|:)=  R22lI[^+("  R^cos0si^+Cos0)09  -  ^cos^sin0^@ 

R-, 

+  (-  T7-sin(t)sinl/r+sin0  )gl 

R2        ^         f  J  (25) 


~  =  R22M  [(sinT/Tcos"(/r-sin20  cosl/rsinYr-  ~cosi/r)  0 
dr  2 

'  2  1  2 

-cos  l/Tsi  -a\fr(p    +(-  <q-  cos0  sinT//'(R-.-R?sinT/r)-cos   -^r 

cos0)  Q0 

R 
-gisin^cosVrQ^rJ  ^^ 


il  =  |_cos(T)sin^  (27) 


Substituting  equations  25,  26  and  27  into  equation  24  and 

:2 


dividing  through  by  the  ( I ,EEl P "~ )  .yields  the  following  ~\jf 


equation  of  motion: 
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R 

(-  -5—  sin0sinl/r+sin0  )Q  +    iff 


+  cos0(l+cos2T^)  Q0 

R 

:2 


+  (^sin2T/r(sin20 -1)+  j—cosl/r)  @2 


+-^sin2T/^  rh  +  #-cosGfosinT/r=  0 

^   R2        y  (28) 


4.   Solving  the  Equations  of  Motion 

The  equations  of  motion  were  solved  for  time  his- 
tory information  of  the  RPV's  position,  velocity  and  accel- 
eration using  the  IBM  Continuous  Si^stem  Modeling  Program 
(CSI-XP)  which  was  available  as  resident  software  on  the  IBM 
360-67  computer  located  in  the  W.R.  Church  Computer  Center, 
Naval  Postgraduate  School.   The  use  of  CSMP  as  a  high-level 
programming  tool  is  quite  similar  to  the  familiar  engineer- 
ing language  of  FORTRAN,  and  programming  details  may  be  found 
in  the  test  by  Speckhart  and  Green,  Ref.  5.   Details  of  the 
actual  program  used  are  presented  in  the  Appendix.   To  use 
this  software  package,  the  three  equations  of  motion  were 

first  solved  for  the  highest  order  derivative  Q  ,  0  and  ty 
in  terms  of  the  three  generalized  coordinates  and  their  first 
time  derivatives.   Then,  in  a  manner  similar  to  an  analog 
computer  solution  of  differential  equations,  the  CSMP  program 

proceeded  to  integrate  these  equations  numerically  for  values 
for  Q  ,  0  ,  "\fr  ,  0  j   0  2.n(3-  ~lp~,    which  were  fed  back  into  the 
equations  for  0  ,  0  and  iff  .   This  iterative  method  of  , 
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solution  was  repeated  at  each  time  increment.   The  resulting 
values  for  0  ,  0  and  ip'    were  listed  as  a  function  of  time 
and  from  these  the  motion  the  RFV  was  studied  to  assess  the 
feasibility  of  the  recovery  system  when  a  particular  set 
of  structural  parameters  and  initial  RPV  conditions  was 
specified. 

Since  the  governing  equations  of  motion  (eqns.  18, 
23  and  23)  were  inertially  cross  coupled,  it  was  necessary 

to  rearrange  the  equations  into  an  inertially  uncoupled  form 
that  would  be  acceptable  for  solution  by  the  CSLEP.   In  matrix 
form,  equations  18,  23  and  28  may  be  expressed  as; 


a 


11 


a 


21 


'31 


L 

where 


12 


0 


32 


"11 


a 


13 


a 


23 


0 


r 


e  -i 


t 


.  0  . 


r  D;o 


=  < 


b 


2  ) 


-  b3  ' 


^  +(_i  _sinT/r)2+cos2T/rsin20 


(29) 


(30a) 


R- 


a12  =  a3l  =  sin0(l-  llsin^r) 


R 


1. 


i^«2. 


L13  =  a21  =  cos0(|pcos^r--g-sin  iff) 


(30b) 
(30c) 


a23  =  cos  y 


a->«  =  1 
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(30d) 
(30e) 


30 


R-.  .  „  R 

.,    =  ~£in0cos'l//-l//^+sin0(^cos'l/r-{rsin2l//')  0 


R 


1 


-{sin  l/rcos  0  -2^±c  o  sty  )Qijf  -sin20cos  1/T00 


R. 


-(l-cos2l/T-2^sinT/r)cos0T^0-  f-^-2  0 


(30f) 


b2   =  -jcos'_T//-sin2  0  0+cos0(l+cos2 1/^)0 0  +sin2?/r'i^0- 


g-sin0cos^T 


R 


b^   =    (-rrsin2T/Tcos^0  -  ^cost/t)  0   -isin2^"0 


cos2T/O60  -  §-cos(pslnty 


(30g) 

■cos0(l+ 

(30h) 


Substituting  equations  30a  through  30h  in  matrix  expression 
29,  6,  0  2-nd  ij/*  were  solved  using  Cramers'  rule,  Ref.  4. 


0  = 


Al 


|  A|     (31a) 


ty  =  112 

1   A 


(31b) 


A. 


0  =  QI 

|  A|     (31c) 


v.'iiere 


A 


a 


11  c'12  c"13 
21  °  a23 
31      a32        ° 


a 


a12   a23   a31-a11   a23+a21   a13 


(32a) 


bl 

a12 

a13 

Al 

= 

b2 

0 

a23 

b, 

1 

0 

-   b3   a12   a23-b1   a23+b2   a^ 


(32b) 


31 


A, 


"11   1   13 


'21   D2  a23 


a31   ^3    ^ 


oil  i  ap "}""  ?al  3  ""  ^  al 1  a"?  ^— ^21^1^ 

(32c) 


A- 


all   a12   bl 


a21 


a31 


0   b. 


1   ~o- 


=  a-j-i  (boa-,  p)-(bPa-)-|-b1api  )+b-j 


-31^u2a12 


(-a 


21a12 


) 


(32d) 


Equations  30a  through  32d  were  programmed  on  the  IBM  360-67 
computer  using  CSI.IP,  the  details  of  which  are  described  in 
the  Appendix.   The  program  evaluated  8  >  0  »  "l/r  ,  Q   i  0  i  ty  j 

0,  0  end' iff    as  a  function  of  time  from  the  initial  con- 
dition of  RPV  impact.   The  evaluation  of  the  computer 
generated  time  history  results  is  described  in  the  next 
section. 
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III.   SYSTEM  DESIGN  ANALYSIS  AND  RESULTS 

A.   UNDAMPED  SYSTEM 

The  reference  or  base  configuration  of  the  recovery 
system  described  in  section  II. A.  and  Fig.  3  was  the  first 
configuration  considered  in  the  analysis.   From  this  con- 
figuration each  parameter  was  varied  individually  for  various 
RPV  approach  speeds  to  assess  the  system1 s  sensitivity  to 
their  changes. 

From  the  outset  of  the  investigation  it  was  felt,  as 
supported  by  early  results,  that  the  motion  of  the  RPV 
along  the  0  coordinate  would  be  a  primary  system  constraint. 
A  guideline  value  for   0 maximum  (  0max)  was  chosen  to  be  90 
degrees,  which  corresponds  to  the  recovery  rod  (R~)  reaching 
the  horizontal.   For  the  recovery  system  to  pass  its  first 
requirement,  an  RPV  approaching  at  50  lenots  would  have  to 
be  recovered  with  the  maximum  amplitude  of  0  not  exceeding 
90  degrees. 

The  results  of  modeling  the  reference  conf iguration 
suggested  that  energy  disipation  would  have  to  be  included 
to  bring  0max  down  to  90  degrees  or  less.   Without  any 
damping  the  system  exhibited  unacceptably  high  magnitudes 
for  0  that  corresponded  to  the  RPV  completely  rotating 

about  the  horizontal  member  when  approach  speeds  exceeded 
25  knots.  This  can  be  seen  by  referring  to  Table  I  where 
0  is  listed  for  the  base  configuration  as  a  function  of 
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time  after  impact  for  approach  speeds  ranging  from  7  knots 
(  0O=2  rad/sec  Ref.  Equation  1)  to  50  knots  (  0O=14  rad/sec). 
Note  that  for  0  =7   rad/sec,  or  approximately  a  25-knot 
approach  speed,  the  magnitude  of  0  at  0.7  seconds  after 
impact  was  greater  than  130  degrees.   The  RPV  has  started 
to  rotate  over  the  horizontal  member.   Variation  of  the 
system  parameters  did  not  correct  this  deficiency. 

3.   DAMPING  ADDED  TO  THE  HORIZONTAL  MEMBER 

In  an  effort  to  correct  the  excessive  RPV  swing  in  the  0 

coordinate  direction,  a  Rayleigh  dissipative  energy  term, 
D,  was  introduced  into  Lagrange's  equation: 

d_  I   3  T  \   3T       1  V      3D_n 

^\^±y   3q±      3q±      >q±    -.  (33) 

•  2 
where   D=4-G-,  0 

This  additional  term  precipitated  through  the  equation 

of  motion  derivation  of  section  II. C.  with  the  addition  of 
only  a  On  9  term  to  the  Q   equation  of  motion,  equation  13. 
This  change  was  easily  incorporated  into  the  computer  pro- 
gram as  indicated  in  the  Appendix. 

An  initial  estimate  of  the  magnitude  of  the  damping 
constant  C-,  was  obtained  using  the  simplified  linearised 
pendulum  equation  below.   The  damping  value  was  varied  to 
detect  the  system's  sensitivity  to  C,  • 

6+2Zoon6  +  wn2Q=  o  (34) 


■5 


where 

'  (S)    2=  £- 
n       Io 

0 

For  values  of  I  =1159  lb-ft-sec2,  K=50  ft-lb  and  £  =  0.15, 
which  was  selected  to  introduce  just  a  small  amount  of 
damping,  the  initial  estimate  of  C-,  was  75  ft-lb-sec. 

The  results  of  adding  damping  C-,  to  the  pivot  of  the 
horizontal  member  in  the  base  configuration  indicated  that 

it  was  only  partially  successful  in  reducing  0max.   This 
can  be  seen  by  referring  to  Table  II  where  0  is  listed  for 
this  configuration  as  a  function  of  time  after  impact  for 
approach  speeds  ranging  from  7  knots,  where  0  =  2  rad/sec, 
to  50  knots,  where  0  =14  rad/sec.   Note  that  the  magnitude 
of  0  was  still  approaching  180  degrees  for  approach  speed 
initial  conditions  in  excess  of  35  knots  (  0  =10  rad/sec). 
Variation  of  system  -oararneters  including  G-.  did  little 
to  further  reduce  0max,  as  shown  in  Figures  4  through  8, 

where  the  base  configuration  of  R-j_=23  ft,  Rp=6  ft,  1/1=  4.66 

2 

slugs,  I  =11^9  ft-lbs-sec  ,  K=50  ft-lb  and  C,=75  ft-lb-sec 

is  plotted  with  only  the  indicated  parameter  varied,  ceteris 
paribus. 

G.   DAMPING  ADDED  TO  RECOVERY  ROD 

In  a  further  effort  to  reduce  the  maximum  swing,  0:nax, 
of  the  recovery  system,  damping  was  also  added  to  the 
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recovery  rod.   The  added  damping  was  included  in  the  system's 
model,  without  cross  coupling,  by  modifying  the  Rayleigh 
dissipative  energy  function  to  the  form  of: 

D=iC162+lrG202  (35) 

The  additional  dissipation  term  precipitated  through  the 
equation  of  motion  derivation  described  in  section  II. C.  as 
the  addition  of  a  C„0  term  to  the  0  equation  of  motion, 
equation  23.   This  change  was  also  easily  incorporated  into 

the  computer  program  as  indicated  in  the  Appendix. 

An  estimate  of  Cp  was  obtained  in  a  manner  similar  to 
that  described  in  the  previous  section  using  the  simple 
pendulum  equation: 

0  +2£o;n0+  (Jn0  =   0  (36) 


where 


2CC.=  °2 


<"n= 


HR22 


vT 


Only  a  slight  amount  of  damping,  C2=3.24  ft-lb-sec, 
added  to  the  recovery  rod  of  the  base  configuration,  in 
addition  to  the  horizontal  member  damping,  was  sufficient 

to  reduce  0max  to  approximately  90  degrees  at  RPV  approach 
speeds  of  50  knots  (  0O=14  rad/sec).   This  can  be  seen  by 
referring  to  Table  III  where  0  is  listed  as  a  function  of 
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time  after  impact  for  approach  speeds  ranging  from  7  knots 
(  0O=2  rad/sec)  to  50  knots  (  0  =14  rad/sec).   The  maximum 

magnitude  of  6  and  Iff   for  the  modified  configuration  at  the 
50-knot  RPV  approach  speed  was  50  degrees  and  63  degrees, 
respectively. 

Further  study  of  the  modified  base  configuration  was 
performed  to  obtain  the  magnitude  of  the  RPV's  linear 
acceleration.   The  absolute  value  of  the  RPV's  linear 
acceleration  was  obtained  using  the  following  equation: 


Acce 


x  +y  +z  (37) 


where  x,  y  and  z  are  defined  by  equations  7a,  b  and  c.   A 
plot  of  equation  37  as  a  function  of  time  for  the  modified 
base  configuration  is  presented  in  Figure  9.   It  may  be 
noted  that  the  magnitude  reduces  to  below  one  "g"  within 
one  second  after  the  RPV  contacts  the  recovery  system. 

The  sensitivity  of  0max  to  each  of  the  system  para- 
meters is  indicated  in  Figures  10  through  16,  where  the 
base  configuration  of  It,  =28  ft,  R?=6  ft,  il=4.66  slugs, 
I0=1159  ft-lb=sec2,  K=50  ft-lb,  C-,=75  ft-lb-sec  and  C2= 
3.24  ft-lb-sec  is  plotted  with  only  the  indicated  parameter 
varied,  ceteris  paribus.   It  can  be  noted  from  Figure  12 
that  the  system  was  sensitive  to  the  recovery  rod  length, 

Rp,  which  nay  be  interpreted  as  an  implication  that  the 
engagement  point  for  RPV  capture  on  the  vertical  cables 
has  considerable  importance  to  system  dynamic  response. 
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The  system  was  also  sensitive  to  the  value  of  C-,  ,  but 
more  sensitive  to  Cp.   This,  however,  should  not  present  a 
problem,  as  C-.  and  C_  will  be  set  by  the  design  and  will 

not  vary  in  operation. 

The  system  was  relatively  insensitive  to  variations  of 

j 

the  point  along  the  horizontal  member,  R-i  ,  where  the  RPV 

was  caught.   The  system  was  also  relatively  insensitive  to 

variations  in  the  moment  of  inertia,  I  ,  of  the  horizontal 

o 

member  and  the  mass,  LI,  of  the  RPV.   The  maximum  displace- 
ment of  0  was  not  affected  at  all  by  changes  in  the  spring 
constant,  K. 

D.   IMPLEMENTING  DAMPING  ON  A  RECOVERY  CABLE 

The  damping  desired  at  the  pivot  of  the  horizontal 
member  can  be  readily  implemented  in  actual  construction. 
However,  it  was  not  readily  apparent  how  the  damping  of 
the  vertical  rod,  which  was  a  cable  in  the  actual  system, 
could  be  implemented. 

Two  simple  pendulum  designs  were  analysed  and  modeled 
using  the  same  techniques  described  to  study  the  recovery 
system.   The  intent  was  to  see  if  pivotal  damping  could 
be  applied  effectively  to  control  the  swing  of  a  cable. 

The  first  design  is  shown  in  Figure  17  and  simulates 
a  cable  design  for  the  vertical  massless  recovery  rod  that 

was  used  in  the  recovery  system  analysis.   With  no  pivotal 
damping,  an  initial  velocity  was  applied  to  M«  to  make 

(£>     approach  a  high  angle  without  exceeding  180  degrees. 

^2max  - J-  °    ° 
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Figure  17   Pendulum  Design  A 
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Then  damping  was  applied  to  the  pivot  of  the  cable  support 
LL.  to  judge  its  effect  on  reducing  the  magnitude  of  0?. 
The  results  are  plotted  in  Figure  18.   This  approach  of 
applying  pivotal  damping  did  not  effectively  reduce  02   v 
as  desired. 

Another  design,  shown  in  Figure  19 »  was  also  analysed 
to  see  if  it  could,  through  pivotal  damping,  control  the 
maximum  magnitude  of  02»   Again  the  design  simulated  a 
cable  instead  of  the  massless  recovery  rod  used  in  the 
analysis.   A  spring  was  used  at  the  pivot  to  give  f/L  a 

near  horizontal  static  deflection.   V/ith  no  damping  applied, 
an  initial  velocity  was  found  for  M2  that  caused  02max  ^° 
approach  180  degrees.   Pivotal  damping  was  then  applied 
to  M-.  to  observe  its  effect  on  the  magnitude  of  02.   The 
results  shown  in  Figure  20  indicate  that  pivotal  damping 
applied  to  this  pendulum  design  was  also  ineffective  in 
reducing  02max. 

While  neither  of  these  designs  demonstrated  effective 
control  of  0?,  they  by  no  means  exhaust  the  possible 
methods  of  implementing  damping  in  a  cable  design,  and 
further  study  is  possible. 
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Figure  19   Pendulum  Design 
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IV.   CONCLUSIONS  AND  RE  CQ;/U,]E  NIDATIONS 

For  the  carousel  recovery  system  to  be  feasible  at 
approach  speeds  of  50  knots,  a  relatively  small  amount  of 
damping  was  required  at  the  pivot  of  the  recovery  rod  in 
addition  to  damping  at  the  pivot  of  the  horizontal  member. 
The  damping  at  both  pivots  was  required  to  prevent  the 
maximum  swing  of  the  point  mass  RPV  from  exceeding  approx- 
imately 90  degrees  in  the  0   direction,  the  coordinate  that 
exhibited  the  largest  change  in  magnitude. 

The  recovery  system  was  sensitive  to  the  vertical 
height  of  the  RPV  at  the  time  it  was  captured.   The  greater 
the  distance,  R~,  between  the  "RPV  and  the  horizontal  member, 
the  greater  the  swing  of  the  RPV  relative  to  the  horizontal 
member.   When  selecting  parameter  values  such  as  the  damp- 
ing constants  under  the  most  severe  conditions,  the  RPV 
should  be  assumed  to  have  contacted  the  recovery  system  at 
the  maximum  length  of  the  recovery  rod  or  cables. 

The  recovery  system  was  relatively  insensitive  to 
variation  in  the  distance  the  RPV  was  from  the  vertical 
mast  when  it  was  caught.   The  system  was  also  relatively 
insensitive  to  variations  in  the  mass  of  the  RPV  and  the 
moment  of  inertia  of  the  horizontal  member. 

Variations  in  the  spring  constant  for  the  torsional 
spring  located  at  the  pivot  of  the  horizontal  member  had 
little  effect  on  the  swing  of  the  RPV.   Since  it  was  not 
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important  to  controlling  the  motion  of  the  RPV,  it  could  be 
removed  from  the  design. 

The  magnitude  of  linear  acceleration  exhibited  by  the 
RPV,  with  the  recovery  system  design  that  incorporated  both 
horizontal  member  and  recovery  rod  damping,  was  within 
Aquila  load  limits. 

Based  on  this  analysis  the  carousel  recovery  system  is 
still  viable;  however,  it  is  apparent  that  the  system  will 
be  sensitive  to  approach  techniques  used  during  the  recovery. 

It  is  recommended  that  this  analysis  be  continued,  in 
order  to  treat  the  RPV  as  a  distributed  mass  with  its  own 
moment  of  inertia  and  coordinate  system,  and  caught  by  two 
recovery  rods  or.  cables,  one  for  each  wing  tip,  instead  of 
one.   Expanding  the  scope  of  the  recovery  system  model  would 
provide  a  more  detailed  analysis  of  the  RPV's  linear  lateral, 
longitudinal  and  normal  accelerations  and  hence  insure  that 
they  were  individually  within  the  RPV's  design  limits. 

It  is  also  recommended  that  the  effects  of  an  RPV 

approaching  at  angles  to  the  recovery  system  be  investi- 
gated; this  should  result  in  increased  motion  in  the  ~\p~ 
direction. 

It  is  further  recommended  that  continued  study  should 
be  made  of  a  means  of  applying  damping  to  the  swing  of  the 
vertical  recovery  cables. 
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APPENDIX 

MODELING  THE  EQUATIONS  OP  MOTION  USING  CSHP 

An  application-oriented  software  package  called  the 
Continuous  System  Modeling  Program  (CSI/IP)  was  used  to  model 
and  solve  the  nonlinear,  cross-coupled  equations  of  motion. 
CSI.IP  is  a  powerful  software  package  that  can  handle  this 
type  of  problem  routinely.   In  the  analysis  performed  in  this 
thesis  two  integrations  using  a  fourth  order  variable  inte- 
gration step  Runge-Kutta  technique  on  each  of  the  three 
equations  of  motion,  Ref.  5.   This  sophisticated  integration 
technique  has  the  advantage  of  automatically  adjusting  the 
time  increment  of  integration  to  meet  the  demands  of  the 
dynamic  conditions  of  the  simulation.   The  absolute  value  of 
the  estimated  integration  error  and  the  relative  magnitude 
of  the  estimated  error  are  compared  with  user  specified 

error-bounds.   The  step  size  is  then  adjusted  to  meet  the 
desired  error  criteria.   In  the  analysis  performed  in  this 

thesis  the  default  value  of  .0001  for  both  the  absolute  and 

relative  errors  was  used. 

The  program  structure  of  CSMP  is  composed  of  three 
segments;  Initial,  Dynamic,  and  Terminal.   Generally,  data 

statements  will  appear  in  the  Initial  segment.   In  addition, 
calculations  that  are  required  to  be  performed  only  one 
time  during  the  simulation  can  be  conveniently  placed  in 
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this  segment.   The  Dynamic  segment  is  usually  composed  of 
structural  statements  that  describe  the  dynamics  of  the 
system  or  explicitly  describe  a  set  of  differential  equa- 
tions.  The  Terminal  segment  is  the  last  segment  in  the 
program  and  is  usually  made  up  of  output  control  statements, 

The  actual  program  used  in  this  analysis  for  the  un- 
damped system  will  now  be  described.   List  of  symbols  used 
in  the  Computer  Programs: 

Rl    -  Length  of  horizontal  member  in  feet 

R2    -  Length  of  recovery  rod  in  feet 

M     -  Mass  of  RPV  in  slugs 

HIO   -  Moment  of  inertia  of  horizontal  member  in  ft-lb- 

2 
sec   about  its  pivot  point 

K     -  Spring  constant,  in  ft-lb/rad, . of  spring  at  pivot 

of  horizontal  member 
X5I   -  Initial  angular  velocitjr  of  the  recovery  rod  in 
the  0  direction  in  units  of  rad/sec  also  equal 
to  approach  velocity  of  RPV  in  ft/sec  divided 

by  R2 

o 
-  Acceleration  of  Gravity  32.2  ft/sec" 


G 

-  Aci 

XI 

-  e 

X2 

-0 

X3 

-t 

X4 

-  6 

X5 

-  0 

X6 

-  t 

X4B0T  -  0 
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X5D0T  -  0 

x6dot  -  i[r 

CO    -  C.  ,  damping  constant  in  ft-lb-sec  for  the  damper 
located  at  the  pivot  of  the  horizontal  member 

002   -  Op,  damping  constant  in  ft-lb-sec  for  the  damper 
located  where  the  recovery  rod  attaches  to  the 
horizontal  member 

AD    -  Determinant  of  A  as  defined  in  Section  II. 0.4 

AD1  -  Determinant  of  Al  as  defined  in  Section  II. G. 4 

AD2  -  Determinant  of  A2  as  defined  in  Section  II. 0.4 

AD3  -  Determinant  of  A3  as  defined  in  Section  II. 0.4 

Bl  -  is  b-,  defined  in  Section  II. G. 4 

B2  -  is  b2  defined  in  Section  II. 0.4 

33  -  is  b->  defined  in  Section  II. C. 4 

Al  -  is  a-,-,  defined  in  Section  1 1. 0.4 


A2 


-  is  a-,  p  defined  in  Section  II. 0.4 


A3    -  is  a-,  ^  defined  in  Section  II. 0.4 

A4    -  is  a0-,  defined  in  Section  II. C. 4 

A,  B,  C,  D,  E,  F,  H,  L,  S,  T  -  were  used  to  simplify 

program  notation  and  are 

defined  in  the  program 
A  description  of  the  program  shown  in  Figure  21,  follows: 

Lines  1-3:     Make  up  the  initial  segment,  parameters 
that  would  remain  constant  during  the 
integration  were  assigned  here. 

Lines  4-34:    Make  up  the  Dynamic  segment. 

Lines  5-13:    Simplifying  notation. 
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Lines  14-28:   The  algebra  used  to  solve  for  0,  0  and  iff 

prior  to  integration  as  described  in 

Section  II. G. 4. 
Lines  29-31:   Perform  the  first  integration  from  angular 

accelerations  to  angular  velocities. 
Lines  32-34:   Perform  the  second  integration  from  angular 

velocities  to  angular  position. 

Lines  35-41:   Make  up  the  Terminal  segment. 

Line  36:      Limits  the  time  period  that  will  be  modeled. 

Lines  38-39:   Prints  output  as  a  function  of  time. 

When  damping  was  added  to  the  horizontal  member  the 
only  changes  that  occurred  to  the  program  were  the  following: 

Line  2:    add   ,00  =  75 

Line  23:   add   -CO  *  X4 - 

When  damping  was  also  added  to  the  recovery  rod  the 
additional  changes  that  occurred  to  the  program  were  as 
follows: 

Line  2:    add   ,002  =  3.24 

Line  24:   add   -002  *  X5 
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